# KRR: compare randomly selected AA and WA

## Preparation

For collection of phenotypes and data split, see `KRR_AAvsWA.md`.

For each site, select all AA and random WA with the same size. If #allAA > #allWA for some site, then keep all WA and randomly select same size of AA. Then assign the selected AA and WA for each data split.

```matlab
ls_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race/scripts/lists';
ABCD_select_allAA_randWA(...
    fullfile(ls_dir, 'phenotypes_pass_rs.txt'), ...
    'race', 'site', fullfile(ls_dir, 'behavior_list.txt'), ...
    fullfile(ls_dir, 'subjects_pass_rs_pass_pheno.txt'), ...
    '/home/jingweil/storage/MyProject/fairAI/ABCD_race/mat/matchANDsplit/20200719', '_pass_rs_pass_pheno')
```

## ------ KRR: regress age, sex, FD, DVARS, ICV, parental education from both behaviors and RSFC ------

### Run KRR; get predictable behaviors (see `KRR_AAvsWA.md`)

### Compute accuracy metric per race group

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_KRR_pCOD_allAA_vs_randWA(...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC', 'lists', 'pCOD_predictable.txt'), ...
        fullfile(proj_dir, 'scripts', 'lists', 'subjects_pass_rs_pass_pheno.txt'), ...
        fullfile(proj_dir, 'mat', 'matchANDsplit', '20200719'), ...
        '_pass_rs_pass_pheno', 120, ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat'))
    ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_KRR_corr_allAA_vs_randWA( ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC', 'lists', 'corr_predictable.txt'), ...
        fullfile(proj_dir, 'scripts', 'lists', 'subjects_pass_rs_pass_pheno.txt'), ...
        fullfile(proj_dir, 'mat', 'matchANDsplit', '20200719'), ...
        '_pass_rs_pass_pheno', 120, ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat') ...
    )
    ```

### Permutation test of accuracy difference

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_PermTest_AAvsWA( ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC', 'lists', 'pCOD_predictable.txt'), ...
        'predictive_COD', ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat') )
    ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_PermTest_AAvsWA( ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC', 'lists', 'corr_predictable.txt'), ...
        'corr', ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat') )
    ```

### Violin plot

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_violin_withnull_AAvsWA(...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC', 'lists', ...
      'pCOD_predictable.txt'),  fullfile(proj_dir, 'models', 'KRR', '20200721', ...
      'reg_AgeSexMtIcvPEduc_fr_y_FC', 'lists', 'pCOD_predictable_colloquial.txt'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat'), ...
      'predictive_COD', ...
      fullfile(proj_dir, 'figures', 'AAvsWA', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC'), 'pCOD_allAA_vs_randWA')
   ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_violin_withnull_AAvsWA(...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC', 'lists', ...
      'corr_predictable.txt'),  fullfile(proj_dir, 'models', 'KRR', '20200721', ...
      'reg_AgeSexMtIcvPEduc_fr_y_FC', 'lists', 'corr_predictable_colloquial.txt'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_FC_allAA_vs_randWA_predictable.mat'), ...
      'corr', ...
      fullfile(proj_dir, 'figures', 'AAvsWA', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y_FC'), 'corr_allAA_vs_randWA')
    ```

## ------ KRR: no confound regression ------

### Run KRR; get predictable behaviors (see `KRR_AAvsWA.md`)

### Compute accuracy metric per race group

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_KRR_pCOD_allAA_vs_randWA(...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'no_reg'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'no_reg', 'lists', 'pCOD_predictable.txt'), ...
        fullfile(proj_dir, 'scripts', 'lists', 'subjects_pass_rs_pass_pheno.txt'), ...
        fullfile(proj_dir, 'mat', 'matchANDsplit', '20200719'), ...
        '_pass_rs_pass_pheno', 120, ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat'))
    ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_KRR_corr_allAA_vs_randWA( ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'no_reg'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'no_reg', 'lists', 'corr_predictable.txt'), ...
        fullfile(proj_dir, 'scripts', 'lists', 'subjects_pass_rs_pass_pheno.txt'), ...
        fullfile(proj_dir, 'mat', 'matchANDsplit', '20200719'), ...
        '_pass_rs_pass_pheno', 120, ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat') ...
    )
    ```

### Permutation test of accuracy difference

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_PermTest_AAvsWA( ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'no_reg', 'lists', 'pCOD_predictable.txt'), ...
        'predictive_COD', ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_pCOD_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat') )
    ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_PermTest_AAvsWA( ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'no_reg', 'lists', 'corr_predictable.txt'), ...
        'corr', ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_corr_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat') )
    ```

### Violin plot

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_violin_withnull_AAvsWA(...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'no_reg', 'lists', ...
      'pCOD_predictable.txt'),  fullfile(proj_dir, 'models', 'KRR', '20200721', ...
      'no_reg', 'lists', 'pCOD_predictable_colloquial.txt'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_pCOD_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat'), ...
      'predictive_COD', ...
      fullfile(proj_dir, 'figures', 'AAvsWA', 'KRR', '20200721', 'no_reg'), 'pCOD_allAA_vs_randWA')
   ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_violin_withnull_AAvsWA(...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'no_reg', 'lists', ...
      'corr_predictable.txt'),  fullfile(proj_dir, 'models', 'KRR', '20200721', ...
      'no_reg', 'lists', 'corr_predictable_colloquial.txt'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_corr_pass_rs_pass_pheno_no_reg_allAA_vs_randWA_predictable.mat'), ...
      'corr', ...
      fullfile(proj_dir, 'figures', 'AAvsWA', 'KRR', '20200721', 'no_reg'), 'corr_allAA_vs_randWA')
    ```

## ------ KRR: regress age, sex, FD, DVARS, ICV, parental education from behaviors ------

### Run KRR; get predictable behaviors (see `KRR_AAvsWA.md`)

### Compute accuracy metric per race group

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_KRR_pCOD_allAA_vs_randWA(...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y', 'lists', 'pCOD_predictable.txt'), ...
        fullfile(proj_dir, 'scripts', 'lists', 'subjects_pass_rs_pass_pheno.txt'), ...
        fullfile(proj_dir, 'mat', 'matchANDsplit', '20200719'), ...
        '_pass_rs_pass_pheno', 120, ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat'))
    ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_KRR_corr_allAA_vs_randWA( ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y', 'lists', 'corr_predictable.txt'), ...
        fullfile(proj_dir, 'scripts', 'lists', 'subjects_pass_rs_pass_pheno.txt'), ...
        fullfile(proj_dir, 'mat', 'matchANDsplit', '20200719'), ...
        '_pass_rs_pass_pheno', 120, ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat') ...
    )
    ```

### Permutation test of accuracy difference

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_PermTest_AAvsWA( ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y', 'lists', 'pCOD_predictable.txt'), ...
        'predictive_COD', ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat') )
    ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_PermTest_AAvsWA( ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y', 'lists', 'corr_predictable.txt'), ...
        'corr', ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat') )
    ```

### Whisker plot

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_whisker_AAvsWA(...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y', 'lists', ...
      'pCOD_predictable.txt'),  fullfile(proj_dir, 'models', 'KRR', '20200721', ...
      'reg_AgeSexMtIcvPEduc_fr_y', 'lists', 'pCOD_predictable_colloquial.txt'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat'), ...
      'predictive_COD', ...
      fullfile(proj_dir, 'figures', 'AAvsWA', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y'), 'pCOD_allAA_vs_randWA')
   ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_whisker_AAvsWA(...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y', 'lists', ...
      'corr_predictable.txt'),  fullfile(proj_dir, 'models', 'KRR', '20200721', ...
      'reg_AgeSexMtIcvPEduc_fr_y', 'lists', 'corr_predictable_colloquial.txt'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_y_allAA_vs_randWA_predictable.mat'), ...
      'corr', ...
      fullfile(proj_dir, 'figures', 'AAvsWA', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_y'), 'corr_allAA_vs_randWA')
    ```

## ------ KRR: regress age, sex, FD, DVARS, ICV, parental education from RSFC ------

### Run KRR; get predictable behaviors (see `KRR_AAvsWA.md`)

### Compute accuracy metric per race group

1. metric predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_KRR_pCOD_allAA_vs_randWA(...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC', 'lists', 'pCOD_predictable.txt'), ...
        fullfile(proj_dir, 'scripts', 'lists', 'subjects_pass_rs_pass_pheno.txt'), ...
        fullfile(proj_dir, 'mat', 'matchANDsplit', '20200719'), ...
        '_pass_rs_pass_pheno', 120, ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat'))
    ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_KRR_corr_allAA_vs_randWA( ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC', 'lists', 'corr_predictable.txt'), ...
        fullfile(proj_dir, 'scripts', 'lists', 'subjects_pass_rs_pass_pheno.txt'), ...
        fullfile(proj_dir, 'mat', 'matchANDsplit', '20200719'), ...
        '_pass_rs_pass_pheno', 120, ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat') ...
    )
    ```

### Permutation test of accuracy difference

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_PermTest_AAvsWA( ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC', 'lists', 'pCOD_predictable.txt'), ...
        'predictive_COD', ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat') )
    ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_PermTest_AAvsWA( ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat'), ...
        fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC', 'lists', 'corr_predictable.txt'), ...
        'corr', ...
        fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat') )
    ```

### Whisker plot

1. metric: predictive COD

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_whisker_AAvsWA(...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC', 'lists', ...
      'pCOD_predictable.txt'),  ...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC', 'lists', 'pCOD_predictable_colloquial.txt'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_pCOD_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat'), ...
      'predictive_COD', ...
      fullfile(proj_dir, 'figures', 'AAvsWA', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC'), 'pCOD_allAA_vs_randWA')
   ```

2. metric: Pearson's correlation

    ```matlab
    proj_dir = '/home/jingweil/storage/MyProject/fairAI/ABCD_race';
    ABCD_whisker_AAvsWA(...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC', 'lists', ...
      'corr_predictable.txt'),  ...
      fullfile(proj_dir, 'models', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC', 'lists', 'corr_predictable_colloquial.txt'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat'), ...
      fullfile(proj_dir, 'mat', 'AAvsWA', 'KRR', 'sig_corr_pass_rs_pass_pheno_reg_AgeSexMtIcvPeduc_fr_FC_allAA_vs_randWA_predictable.mat'), ...
      'corr', ...
      fullfile(proj_dir, 'figures', 'AAvsWA', 'KRR', '20200721', 'reg_AgeSexMtIcvPEduc_fr_FC'), 'corr_allAA_vs_randWA')
    ```
